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We study the classical stability of an anisotropic space-time seeded by a spacelike, fixed norm, 
dynamical vector field in a vacuum-energy-dominated inflationary era. It serves as a model for 
breaking isotropy during the inflationary era. We find that, for a range of parameters, the linear 
differential equations for small perturbations about the background do not have a growing mode. 
We also examine the energy of fiuctuations about this background in flat-space. If the kinetic terms 
for the vector field do not take the form of a field strength tensor squared then there is a negative 
energy mode and the background is unstable. For the case where the kinetic term is of the form 
of a field strength tensor squared we show that perturbations about the background have positive 
energy at lowest order. 



I. INTRODUCTION 

Inflation has become the standard paradigm for very 
early Universe cosmology. It supposes an era when the 
Universe was dominated by vacuum energy. During this 
era any classical inhomogeneities are smoothed out by 
the exponential expansion of the Universe. The inhomo- 
geneities that we observe today originate in small quan- 
tum fluctuations that have their wavelength rcdshiftcd 
outside the horizon by the exponential expansion of the 
Universe. At some time the Universe exits the inflation- 
ary era and transitions to a radiation dominated era and 
then eventually to a matter dominated era. After in- 
flation, the Universe expands slower than the speed at 
which light travels. This allows the perturbations with 
wavelengths that are outside of the horizon to eventually 
reenter the horizon and generate the inhomogeneities we 
observe in the microwave background (as well as the large 
scale structure of the Universe) . 

Since we have no direct probes of the inflationary era, 
we are compelled to consider the possibility that some 
tenets of physics — tenets that are fundamental to our 
current understanding of the Universe — only reflect post- 
inflationary developments. One such tenet is that rota- 
tional invariance is not spontaneously (or explicitly) bro- 
ken. This symmetry implies angular momentum conser- 
vation and is necessary for the classification of elementary 
particles by their spin (i.e., their angular momentum in 
their rest frame). 

Recently the possibility that rotational invariance is 
broken during the inflationary era has been studied [1- 
6]. In Ref. [1] it was noted that if the breaking of ro- 
tational invariance is small there is a very simple and 



unique signature on the spectrum of density perturba- 
tions and hence on the anisotropy of the microwave back- 
ground radiation. Ackerman et. al. also introduced a 
simple model for breaking rotational invariance during 
the inflationary era. They calculated the influence of 
the breaking of rotational invariance on the spectrum of 
density perturbations, verifying the expectations based 
on their general arguments. The Lagrange density for 
the model has the usual Einstein-Hilbert action, a cos- 
mological constant term with vacuum energy density pA 
and the following terms containing the four-vector field 

- /?3V^u"V,u^ + \{u^'u^ ~ m2) . (1) 

Here A is a Lagrange multiplier that enforces the con- 
straint Qf^^u^u'^ — m? . We take rri^ > so the four- 
vector u'^ is spacelike. 

This model has a homogeneous but anisotropic back- 
ground solution to the classical equations of motion. We 
choose the xa-axis to be aligned along the four-vector 
field, 

u° = 0, u^=u^ = Q and ^ (2) 

b(ty ^ ' 

resulting in a space-time metric of the form, 

ds^ = -At^ + a(<)2dxi + h{tfAx\. (3) 

The breaking of rotational invariance by the four-vector 
field causes the a;3-direction to expand at a different rate 
than the x\ and X2 directions. Explicitly,^ 

a(i) = e«"*, 6(t)=e^^*, (4) 
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^ It is assumed that the dynamics immediately preceding the in- 
flationary era was rotationally invariant so that a(0) = b(0). 
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where the Hubble parameters are related to the vacuum 

energy, the scale of rotational invariance breaking, and 
parameters in the vector Lagrangian by the following 
equations, 



Ha = -= Hb{l + WwGPim^), 
a 



''~ b~ y {1 + 87rG/3iTO2)(3 + 327rG/3irn2) 



(5) 



Notice that the background solution does not depend on 
the parameters /?2,3 and as m ^ the rate of expansion 
is the same for all spatial dimensions. 

The purpose of this paper is to study the classical sta- 
bility of this solution. We first consider the linear differ- 
ential equations that result from expanding the classical 
equations of motion to linear order about the background 
solution discussed above. We find that, for small G/SjW^, 
small perturbations do not grow provided that /3i > 
and /?i -I- /?2 + /Js > 0. Next we consider the energy of 
these solutions in flat-space. Negative energy modes indi- 
cate an instability that might not show up in the analysis 
of the linearized equations of motion. 

A special role is played by the case Pi + + Ps = 
where, in fiat-space, the kinetic terms for the vector field 
take the form of a field strength tensor squared. In this 
case, we find that when the energy density is evaluated 
to quadratic order in fluctuations about the background, 
all propagating modes have positive energy. However, 
when that equality is not satisfied, there is a mode that 
propagates with negative energy. 

Models with vector fields that spontaneously break 
Lorentz symmetry were first introduced in [7]. Similar 
fixed-norm timelike Lorentz-violating vector field mod- 
els have been extensively studied; constraints required 
by theoretical and observational consistency have been 
placed on parameters in this theory [7-19]. For a review, 
sec [20]. 

Others have classified gauge invariant perturbations in 
anisotropic scenarios [21, 22]; however, we choose to work 
in a particular (nonstandard) gauge for calculational con- 
venience. 



II. STRATEGY 

The effect of an isotropy-brcaking vector field on the 
expansion of the Universe during inflation was recently 
studied [1]. In this model, the four- vector is non- 
zero only during the time interval < f < f*, where 
t = is the beginning of inflation and is the end of 
inflation. We assume that the dynamics is rotationally 
invariant during reheating and thereafter. During the 
time interval < t < f*, the dynamics of interest in this 



paper is governed by the action 

1 



S = / d'^xs/^ 



ISttG 



R- pK + Cu 



(6) 



where £„ is given in Eq. (1). The homogeneous 
background inflationary space-time solution is given in 
Eqs. (3), (4) and (5). 

The energy- momentum tensor for u** derived from (1) 
is [9], 

T^«) = 2/3i(V^w^V.w^- Vw^VpU.) 

-2[Vp(n(^ J^)) + V,K J(^.)) - V^(u(^J.)0] 
+2m-'^u„VpJ'"'u^u^ + g^uCu, (7) 

where J''o- is the current tensor. 

Given Eqs. (2) and (7), the nonvanishing components 
of the background stress tensor are. 



p(m) 
'33 



(3im^ 



266 - 4 

a 



(8) 



The components of the energy-momentum tensor in our 
chosen background are independent of P2 and (3s. For a 
non-trivial vector contribution to the stress-tensor (and 
consequently for an anisotropic metric resulting from the 
fleld equations), one must require that Pi ^ 0. Note that 
the background solution does not satisfy the weak energy 
condition. 

An inflating background solution requires. 



SnGm^Pi > --. 



(9) 



This is not a strong boimd on Pi since we are interested 
in backgrounds that have a small violation of isotropy, 
i.e., Gni^ « 1. 

In this paper we study the classical stability of the ho- 
mogeneous background solution to the equations of mo- 
tion. We expand the equations of motion to linear or- 
der about the homogeneous background. Since the back- 
ground is homogeneous it is convenient to Fourier trans- 
form the small fluctuations about the background in the 
co-moving spatial coordinates {xi,X2,Xs) and examine 
the time dependence that the linearized equations of mo- 
tion imply. The physical wave vectors fcj , i = 1, 2, and fcs 
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are related to the co- moving wave vectors by fcj = ki/a{t) 
and fca = k^/h{t). 

We begin with a flat-space stabihty analysis neglecting 
gravity since, when H, such an analysis captures the 
essential physics. Then we perform a stability analysis 
including gravity in two regions, k ^ H and k <^ H. 
Wc do not treat the transition region where the physical 
wavelengths of the modes cross the horizon. In the final 
section of this paper we examine the flat-space energy 
of solutions to the equations of motion. Negative energy 
solutions indicate an instability that might not be evident 
from a study of the linearized equations of motion. 

III. STABILITY ANALYSIS NEGLECTING 
GRAVITY 

In the short-wavelength limit, fc ff, it is physically 
intuitive that modes will not be able to resolve space- 
time curvature. Therefore we should be able to derive 
dispersion relations for small fluctuations in it^ about a 
flat background space-time that will receive corrections 
suppressed by Gm? when gravity is included. In this 
case, we have the action, 

S = Jd^x{-(ii{d>'un{d^u^)-02{d^u^'f 

- p3id^u^){d^u^) + \{u^u^-m^)). (10) 

The Euler-Lagrange equation for A gives, 

r?^,M^M" = m^ (11) 

while the equation of motion for the vector field gives, 

Xu" + pid^d'^u'' + {(32 + (33)8'' {d^u") = 0. (12) 

We can use Eq. (11) to solve for A by taking the inner 
product of the vector Euler-Lagrange equation (12) with 
the vector, 

^ = - (^^) ^^''(S.^'^). (13) 

We consider small perturbations about the back- 
ground, ViP" = mrf^, A = 0, and = iq^uj. We denote 
small perturbations about the background vector field by 
u** so that. 

Since we are neglecting gravity, all greek indices are 
raised and lowered by the Minkowski metric. The per- 
turbation about the background value of the Lagrange 
multiplier, (5A = A — A, is a function of the vector field 
perturbations through Eq. (13), 

mJA + (/J2 + /33)%(a.t;") = 0. (14) 



Prom Eq. (11) we see that U3 = 0. And by expand- 
ing Eq. (12) to first order in perturbations wc have the 
following equations for the non- trivial vector field pertur- 
bations (i e {1,2}): 

/3ia,5'^i;° + (/32 + /33)a°(5,t.'') = 0, 
/3i9,aV + (/32 + /33)a'(a,t;'^) = 0. (15) 

We Fourier transform the components of , and in 
terms of the Fourier modes (transformed in both space 
and time), 

?;°(w, fc) = k39i{ui, k), 
and v^{u, k) = ki02{oj, k) + eijkjil){u, k), (16) 

where i,j G {1,2} and wc have arranged for all scalar 
{9) and pseudoscalar^ {tp) components components to 
have the same mass dimension. With this decomposition, 
linear combinations of the Fourier transform of equations 

(15) become, 

[Pi{kl + kl) - (/3i +132+ /33)^'] fca^i 

+ {f32+f33)ukle2 = 0, (17) 

[(3i{kl + kl-u^) + {P2+Pz)kl] 02 

- (/32 + /33)wfc3^1 = 0, (18) 

/3ifci(fci+A;2-a;2^V = 0, (19) 

where fc^ = Yh^i hki. 

We first consider the case where (3i + fi2 + (33 7^ 0. 
The scalar mode 6\ = {k\/ijjk3)d2 propagates with the 
dispersion relation, 

tJ' = k'i + kl. (20) 

The pseudoscalar mode also propagates with this disper- 
sion relation. The scalar mode 9\ = {u/k3)92 propagates 
with the dispersion relation. 

When (3i+ (32+ Pa — 0, Eq. (17) becomes a constraint 
equation and therefore the number of scalar modes is 
decreased by one while the pseudoscalar equation re- 
mains dynamical. A mode disappears at first order when 



^ We use the name "pseudoscalar" for the tp's not because these 
objects are odd under parity but in order to distinguish them 
from the 0's. The two different kinds of scalers decouple. 
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/3i +/32 +/33 = because there is an enhanced gauge sym- 
metry. At first order in perturbations, the transformation 
Vfj, + dfj_a{t, x\ , X2) loaves the Lagrangian invariant. 

The scalar function a is restricted to have no X3 depen- 
dence at first order by the constraint u^u'' = m^. The 
disappearing mode, 

v°{uj, k) = uje2{uj, k), v'{uj, k) = k'02{uj, k), (22) 

satisfies (in coordinate space) v^(t, x) = dn(j){t, x) for ji = 
0, 1, 2. When (32 + (3^ = Eqs. (15) imply that 4>{t, x) 
is independent of .X3. Thus the mode (22) becomes a 
gauge artifact when /3i + /32 + /^a = 0. 



IV. STABILITY ANALYSIS INCLUDING 
GRAVITY 

In this analysis we have four vector field fluctuations 
and ten metric fluctuations to consider. However, we can 
use diffeomorphism invariance to remove some of them. 
Under an infinitesimal coordinate transformation, 
xi^ + the vector field transforms to + Au^ where, 

Au^ = ui'd^^^ - ^"a^u^, (23) 

while the metric transforms to g^^ + ^g^v where, 

^g^u = -(V^e^ + V^^u). (24) 

We choose; a gauge where the first order fluctuations 
in the (contravariant) four-vector field about their back- 
ground, = r}^^m/b{t), vanish {i.e. 5{u'^) = = 
where u'' = + v^^). This gauge condition is satisfied 
by fixing, 

-V^ = Aw^lptordor 



m 

m 



m 

wr 



It is important to note that the gauge freedom has not 
been completely exploited. The condition 5{u'^) = is 
invariant under gauge transformations of the form xi^ — > 
xt^ -\- C^f^ where, 

C = (/o(i,-Tl,.T2),/l(i, Xi,X2),f2{t.,Xi,X2), 

fsit. Xi,X2) - HtX3fo{t,Xi,X2)). (25) 

Once the above fom- fmictions of three coordinate vari- 
ables are chosen, the gaiige has been completely specified. 
This residual gauge freedom will be used in the analysis 
of long-wavelength fluctuations. 

The vector field satisfies the Lagrange multiplier equa- 
tion, = m^. This equation and our choice of gauge 
imply that, 

2 

Til 

{6g^.)uV + 2g^,{vnu^ = (Sgss)-^ = 0, (26) 



and thus we are left with nine metric fluctuations to con- 
sider. 

Without loss of generality, we have chosen coordinates 
such that the background vector field lies entirely along 
the x^ direction. Since isotropy is broken, only an SO (2) 
spatial symmetry remains. Let the Roman indices i,j,k 
be 80(2) indices that run from 1 to 2. Define the 'co- 
moving' metric perturbations, h^^, through the following 
equation, 

ds'^ = -{1 + hoo)df +2a{t) hoidtdx'' 

+ 2h{t)ho3dtdx^ + 2a{t)b{t)hizdx'dx^ (27) 
+ a{tf {5ij + hij)dx'dx^ + b{tf{l + hssXdx^f 

We work to first order in the small perturbations. After 

making our gauge choice, there are nine perturbations, as 
well as nine independent equations for the system. One 
can show that the equations of motion for the vector 
field are equivalent to the energy-momentum conserva- 
tion equations. Furthermore, energy-momentum conser- 
vation (and thus the equations of motion) follow from 
Einstein's equations. The identity to first order in per- 
turbations is, 



where, 



Finally, only nine of Einstein's Equations are indepen- 
dent. The interdependence among Einstein's equations 
in our gauge can be found by working out the identity 
above for u = 3, since the left hand side vanishes identi- 
cally, 



id, + 2iH, + H,))^-^=S^^.^^^''^' 



b{t) a{tYh{t) b{tf 

Thus we take as our complete set of equations at first 
order in small perturbations all of Einstein's equations 

except for 5Et,z- 

We Fourier transform, 

Ku{t,kuk2,k3) = j j^-^hi,4t,x^,x'^,x^)e-'''-^, 
and define, 

ki = ki/a{t), k]_ = kiki, BM.dks = kz/b{t), 

where i G {1,2} and summation is implied over repeated 
indices. Also define the constants, 

A„ = 87rGTO^/3„ and A = SttGpa- (28) 

Then the Fourier-transformed field equations become: 
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^£^00 = ((A + fcgCAi + A2 + A3)) hoo + 2ik3{XiHb - Ha)ho3 - i{Ha + Hi,)kiho, - k^hh^^ 

+ ^ (^i + - 2A2) + (i?a + Hi,)dt) hu - hikhij = 0, (29) 



= -iki{Ha + Hb)hoo - k3kiho3 + (2A - 2HI - 2HaHb - 2H^{1 + Ai) + kl + kl{l + 2Ai + 2A3))/ioi 

a[t) 

- kikjhoj + ik3{-Ha + Hb{l + 4Ai) + dt)hi3 + i{kjdtlHj - kidthjj) = 0, (30) 



^ = ihiXiHb - Ha)hoo + (a - 3Hl + Xiim^Hb + Hi + 2kl) + ^(1 + 2Ai - 2A3)fci) /i03 

- ^hhhoi + iki^{l + 2Ai - 2A3) (i/a -Hb + dt)h^3 - ih^iHa -Hb{l + 2Ai) + dt)h^^ = 0, (31) 

^ = 5ij (^-{Hl + H^Hb + (1 + M)Hl) + ifci + 1(1 - 2A2)fc3' - \{Ha + Hb)d^ hoo 

— —kikjhoo + ik3dij{Ha + (1 + 2Xi)Hb + dt)ho3 — ik(^i{Ha + Hb + i9t)/ij)o + i^ij kl{Ha + Hb + dt)hoi 
+ k3{Sijkihi3 - k^ihj)3) - ^^^(^3(1 - 2A2) + {2H^ + Hb + dt)dt)hkk 

+ {k- {Hi + HaHb + Hlil + Ai)) + + 2Ai + 2A3) + ^{2Ha + Hb + dt)dt^ h^j = 0, (32) 

25 E r, _ _ _ 

= -k3kihoo - iki{l + 2Ai - 2A3)(2if6 + dt)ho3 - ik3{3Ha - Hb{l + 4Ai) + dt)hoi 



a{t)b{t) 

+ (2A + Hl{2Xi - 2X3 - 5) + H^Hb{l + lOAi - 2A3) - 2^^2(1 + Ai - 2A3) + 4^|Ai) hia 

+ (1 + 2Ai - 2A3) {5ij ikl + {2Ha + Hb + dt)dt) - hh) hjs + h {hhjj - kjhtj ) = 0. (33) 



A. Short-wavelength limit of the field equations 



Here we consider the limit (H/k) 0, which corre- 
sponds to modes with physical wavelengths much shorter 
than the Hubble radius. Such modes have periods much 
shorter than the Hubble time; therefore we can treat the 
ki as constants independent of time {i.e. Ha.b arc effec- 
tively zero in this limit). The background solution in 
Eq. (5) gives VX Hb and thus we also take (A/^^) 0. 



The various degrees of freedom in the field equations 
decouple into six equations governing the scalar modes 
and three equations governing the pseudoscalar modes if 



one uses the decomposition {i,j = 1, 2): 



^00 






hoi 


= ^3 (fci02 + e 


ijkj^i) , 


ho3 


= ^303, 




hij 


= ^y^04 + 


k"i k j 3 








hi3 


= fcs (fcjOe -1- e 


^Jk,^3) ■ 



Note that all of the scalar (8) and pseudoscalar (^) com- 
ponents have the same mass dimension. We make the 
ansatz for the time dependence of the fields above, 

e;(i, k) = @i{uj, k)e^^ and ^i{t, k) = ^/(w, k)e''^\ 
The decomposition (34) aids in solving the field equa- 
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tions in the usual way: linear combinations of Einstein's 

equations that resemble the above decomposition lead to 
equations coupling only (pseudo) scalar fields. 



There are six scalar fields in the mode decomposition 
(34). The functional and algebraic constraints on these 
fields are given by the field equations: 



J 



SE, 



00 



SEo3 



= fc|(Ai + A2 + A3)ei + 



^ + f(l-2A2) 



G4 - ^65 - fcfee, 



+ 2Ai + 2A3)e2 



2 



k^ 
4 

2 



kiQf 



h(t)kj 2 



A-2 



^3 
5Eu 
a(t)2fc2 



fc2 

^(l + 2Ai-2A3) + 2fc2Ai 



Lo I ki^ ki 



2A3)ef 



^ + -klil-2X,) 



ei - ^3^62 - 2^63 + 



+ fc|(Ai + 2A2 + A3 - 



64 + klQe, 



kikjSEij fc| ^ ^f^3 rS 

-_^ = -(l-2A2)Oi-^03 



4 4 ^ 



2Ai - 4A2 - 2A3) 



rE2 



64 



^(l + 2Ai+2A3)- — 



LSE, 



i3 



a{t)b{t)k3kl 



2 



2 



62 - ^©1 + 



^^(1 + 2Ai - 2A3)e3 + ^(©4 - ©5) - 



(l + 2Ai-2A3)-2Aifc| 



©5, 

©6. 



Let us consider the dynamical degrees of freedom. 
We shall schematically denote LvQ{u),k) ~ dtQ{t,k) ~ 
Q{t). The first equation is a constraint. When /3i + 
f32 + P3 7^ 0,'^ it can be solved to give &i{t) = 
©i(©4(i), ©5(t), ©6(i)). The next two equations can 
be solved to give ©2(t) = ©2(©4(i), ©5(0> ©6(i)) and 
eait) = ©3(©4(t),©5(t),©6(t))- In the last three equa- 
tions, 02 and 03 appear only as 62 (t) and 63 (t). Thus, 
after making the substitutions for ©1, ©2, and ©3, the 
last three equations form a system of second order or- 
dinary differential equations of three functions. Thus, 
when /3i + /?2 + /J3 7^ 0, we expect three distinct dynam- 
ical modes. 

When /3i + /?2 + /?3 = 0, the function ©1 drops out 
of the constraint equation. Then the first equation 
can be solved for ©4(©5(t), ©6(t)) and the next two for 

02(t) = 02(05W,©6(i)) and ©3(t) = ©3(©5(i),06(i)). 
Since 0i enters the above equations with no time deriva- 
tives, one of the last three equations can be used to elim- 
inate ©1, and the remaining two become second order 
differential equations in time of two functions, ©5(t) and 
©6(t). Thus, when /3i -|- /32 + /33 = 0, we expect only two 
distinct dynamical modes in the scalar sector. 

Here we present the propagating modes and their dis- 
persion relations. The first scalar mode relates the am- 



^ Recall the definition, An s SnGmPfin- 



plitudes of the various fields by, 



oj 2a)^ — kj 



©3 = = ©6, (35) 



and propagates with dispersion relation, 

Lo^ ^kl+kl{l + 2X1 + 2X3). (36) 

The next scalar mode is characterized by the amplitude 
relationships, 



©1 = ©4 = ©5 

2^3^ 



-©3 



2k3Uj 

UJ^+kl 

2 



02 



©6, (37) 



(38) 



fti(l + 2Ai+2A3)"" I + 2A1+2A3 
and propagates with dispersion relation, 

p / Ai(l + 2Ai+2A3; 

The last propagating scalar mode only occurs when 
/3i + /32 + /33 7^ and is characterized by the amplitude 
relationship, 

_ _ 2fc3(l + Ai) 

- IJ^' - o.(l + 2Ai+2A3)°' 



2(1 + Ai) 
1 + 2Ai + 2A3 



-©fi = 



kl 



-©1 



l+Ai 



(kl 



^ fcgAi(l+2Ai+2A3) 
^ 1+Ai 



-04, 



(39) 
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and dispersion relation, 



Ai(l + 2Ai +2A3)\ /I-A1-3A2-A; 



(1 + Ai 



Ai + A2 + A3 



(40) 



This last mode corresponds to the flat-space mode (22). 
As we found in the flat-space analysis, this mode is absent 
when 01+132+ f3s = 0. (Note that in the short-wavelength 
limit, a(t) and b{t) are treated as constants, so covariant 
derivatives become ordinary derivatives and thus /3i -I- 
/32 + /33 = corresponds to the case where the Lagrange 
density takes the form of a field strength squared — as in 
the flat-space analysis.) 

For general Fourier components, we can guarantee that 
^ if and only if the coefficient of the term is 
positive semideflnite. Thus we see that the stability of 
the modes gives the following conditions on the A„'s (and 
thus, recalling Eq. (28), the /J^'s): 



A1 + A3 > 



\2 
-^3 



Ai(l-H2Ai-H2A3 
(1 + Ai) 



Ai(l + 2Ai + 2A3) 
Ai + A2 

1 - 2A2 
Ai -l- A2 -l- A3 



^ 0, 



(41) 



- 1 > 0. 



We see that only the last equation gives us a condition on 
02- In the limit that Gm^ « 1, we have the following 
condition, 



where A„ = 87rGm^/3„. Here, the first equation is a 

constraint and there arc two distinct dynamical modes. 
The pseudoscalar eigenmode {2ujk3/k^)^i = '^2, ^'3 = 
propagates with the dispersion relation. 



io' = ki 



kl{l + 2Xi+2X3). 



(43) 



The second pseudoscalar eigenmode is given by 
(2fc3/w)(l + 2Ai + 2A3)*i = (1 + 2Ai + 2A3)§2 = 2*3 
and propagates with the dispersion relation. 



Ai(l + 2Ai +2A3) 
Ai+A? 



X2 
^3 



(44) 



Note that the value of 02 does not affect the stability of 
any of the propagating pseudoscalar modes. Also, the 

two pseudoscalar dispersion relations were already found 
to characterize two of the propagating modes in the scalar 
sector. 

When /?i = — /33 and 02 = (the Maxwell case), all 
modes propagate with the speed of light. However, all 
modes would not necessarily propagate with the speed of 
light if we only required (3i + 02 + Pa = 0- 

In the Hmit Gm^ << 1 (or equivalently vn? « M^, 
where Mp is the Planck mass) we have A„ << 1. The flat- 
space analysis considered earlier is insensitive to 0(A„) 
corrections to the dispersion relations and (of course) ne- 
glects gravitational degrees of freedom. Here we show 
that the more general analysis simplifies to the flat- 
space analysis. Since, in our gauge, 5(u^) = Sqi^vU" = 
(m/b) Sg^s, comparing Eq. (16) and Eq. (34) we expect 
that, 



Ai 



/3i 



Ai + A2 + A3 



/3l + /32 + 03 



>0, 



(42) 



which is identical to what we found in the flat-space anal- 
ysis. The other two constraints are trivially satisfied in 
this limit. 

The following equations for the pseudoscalar fields fol- 
low from the field equations: 



= 



2eijdEoikj 
kskj^a{t) 

h V 2 

eijkikiSEij 



[fci + fc|(l + 2Ai + 2A3)] *i 



i*2 + ^32^3 



kla{t) 
,2 



2 ^ 2 



eijSEi^k 



fcg 

4 



(H-2A1 +2A3 



*2, 



k3k\a{t)b{t) 



-^l =!^*2 

2 4 



(l-F2Ai-2A3)-2AiA;3 



3, 



6*2 OC 96 



and V' OC 



(45) 



In the A„ << 1 limit, the dispersion relation for the 
scalar mode, (38), and the pseudoscalar mode, (44), has 
the first order approximation. 



ki 



1 



(A1+A3)- 

Ai 



(46) 



with ampHtude relationship (ksUj/kY)^^ = ©e- The 
dispersion relation in Eq. (40) has the first order approx- 
imation. 



— kj_ + A13 



/3i 



Pi +02+ 03 



[1 - 3A2 + A3] , (47) 



with amplitude relationship (k^/uj^ 6)3 = Og. 

The A„ Hmit of dispersion relations (38), (40), 
and (44) and the corresponding mode amplitude rela- 
tionships, along with Eq. (45), lead to the propagating 
modes and dispersion relations foimd in section HI. (In 
the flat-space analysis, one can rescale the Fourier modes 
in order to put them in the form of k± and ^3 since 
this simply results in an overall constant rescaling of the 
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integration measure.) We found only three (two when 

Pi + P2 + = 0) distinct modes in the analysis neglect- 
ing gravity because the = = 83 and vJ'a = modes 
are purely gravitational. The purely gravitational modes 
have the following exact dispersion relation, 



cj"^ = k'i + kl (1 + 2Ai + 2A3) . 



(48) 



As A„ ^ 0, the gravitational modes have the usual gravi- 
ton dispersion relation (as expected). We see that the 
intuitive flat-space analysis is reproduced by the more 
general analysis involving coupling to gravity. 

In [10], Lim carried out a similar perturbative analysis 
of a fixed-norm, timelike, Lorentz-violating vector field 
in a de-Sitter background. Lim considers the model of 
Ref. [9] which is a slight simplification of the model in 
Rcf. [17]."' He found rcscaled mode propagation speeds 
very similar to those that we found when /3i-|-/32+/?3 7^ 0. 
There is a one-to-one correspondence of the inverse of 
his propagation speeds with our aJs-direction propagation 
speeds when —ra^. 



B. Long-wavelength limit of the field equations 

We now consider the behavior of the modes after they 
cross the Hubble Horizon and work in the limit where 
k/H <C 1. First simplifying Einstein's equations in 
this limit and then performing a mode decomposition 
yields simple differential equations for the various modes. 
The field equations decouple into six equations governing 
the scalar (Q) modes and three equations governing the 
pseudoscalar {^) modes if one uses the decomposition 
(i,i = l,2): 



^00 
hoi 
ho3 
hij = 

hi3 = 



ei, 
1 

kl 
63, 

s. 

2 
1 



ikie2 + e^JkJ'^l)^ 



(49) 



^64 



2 



65 + 



k'i 



{kiOe + ^ijkj^s) , 



where k± = ^/kf + fc^- Note that all of the scalar (G) 
and pseudoscalar (*) components have the same mass 
dimension and are zeroth order in k±,k3. 



The scalar mode Einstein equations are: 

= iJf,(3 + 4Ai)ei(t) + ate4(t), 

= 4Hb{l + Xi){Ht,{3 + AXi) + dt)ei{t) 

+ {{3 + AX,)Hi,dt+d^)e4{t), 

= {Hi,{3 + 4Xi)dt+d^)Q5it), 

= {2H^Xi{3 + 2Xi) + Hi,{3 + 4Xi)dt + d^)Qe{t), 

where we have used the background solution, Ha = 
Hb{l + 2Ai) and A = H^{1 + Ai)(3 + 4Ai). The other 
two scalar equations are trivially zero. This means that 
02 (t) and 83(f) are imconstrained fimctions of time. The 
solution for the remaining scalar modes is: 

81(0 = ei(0)e-*^''(3+4Ai)^ (50) 

Qiit) = 0i(O) (e-*^''(3+4Ai) _ ^ Q^^o)^ (5;^) 

85(4) = 85(0)e-*^''(3+4Ai) +^g(e-fH.(3+4Ai) _ 

(52) 

86(0 = e-'*^^^i (86(0)e-3*^' + Aeie-"^*"' - 1)) , 

(53) 

where A4 and are arbitrary constants. One should 
note that Ha-Hb = 2XiHb and i7b(3 + 4Ai) = 2Ha + Hb. 

The nontrivial psuedoscalar mode Einstein equations 
imply, 

= {Ht,{3 + 4Xi)dt+d^)^2it), 

= {2H^Xi{3 + 2Xi)+Hb{3 + 4Xi)di + df)^3it), 

where we have again used the background solution, 
Eq. (27). The third equation vanishes identically, so there 
are no constraints on ^'i(t). 

The solutions for the other modes are, 

*2(t) = *2(0)e-*^^(3+4AO ^ 5^ j^g-t//,(3+4Ai) _ ^ 

(54) 

^^/3{t) = e-2*«''^i (*3(0)e-3*^^ + 53(6-3*^" - 1)) , 

(55) 

where B2 and S3 are arbitrary constants. 

The unconstrained scalar and pseudoscalar modes, 
82 (t), 83 (t), and 4'i(i), are, unsurprisingly, gauge arti- 
facts. Importantly, ^3{t) is unchanged under any trans- 
formation of the form (25). We show this explicitly in 
the appendix. 

That ^'3(t) is not a gauge artifact implies that /3i (re- 
calling that Ai = SnGm^lSi) must be greater than zero. 



/?! >0, 



(56) 



The action in [17] includes a term quartic in the vector field. 



in order for this mode to decay. Thus, assuming Eq. (56), 
no comoving modes grow in the limit k/H << 1. 
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For the timelike vector [10], the /3i > bound guaran- 
tees a positive definite flat-space Hamiltonian at second 
order in perturbations; however, /3i > is not a require- 
ment for lowest order perturbative stability. For the case 
of a timclike vector in de-Sitter space, all comoving modes 
decay as l/a{t) [10]; however, for a spacelike vector we 
find that &e{t) and ^'3(t) fall off more slowly than l/a{t) 
while 6i(f), 64 (t), 65 (t) and 'i>2{t) decay more quickly. 
(Here we assume that Ha — = 2iJ{,Ai is smaller than 

Ha.) 



V. FLAT-SPACE ENERGIES 

In the previous sections we have studied the classical 
equations of motion linearized about the homogeneous 
background solution. For small Gm? they were stable 
provided /?i > and /Ji -|- /32 + /?3 > 0. In this section we 
examine the classical energy of solutions to the flat-space 
equations of motion. Negative energy solutions indicate 
an instability that might not have shown up in our study 
of the linearized equations of motion. 

In flat-space one can write the vector fleld action as, 

Su = - j (fxl^^F^^F^'^ + (/3i + /32 + /33) {d^u^ f 

-Xiu^u^" -w?)), (57) 

where = d^Ui, — di,u^ and boundary terms are set to 
zero. Let i, j G {1, 2} as in previous sections. The energy 
functional derived from the above action is,^ 

- y d^x\{u^u'^ - m^), (58) 

where, in the last expression, all indices are intentionally 
lower and repeated indices should be summed without 
factors of the metric. 

We take = {vQ,vi,V2,m -\- v^) as in section III. 
Physical field configurations must satisfy the Lagrange 



^ The Lagrange multiplier constraint is holonomic. Whether a 
holonomic constraint is imposed before or after forming the 
stress-energy tensor has no effect on the tensor's final form. 



multiplier equation of motion, = m^, which implies, 

Thus V3 is actually second order in perturbations. This 
expansion also makes manifest that an expansion order- 
by-order in perturbations is equivalent to an expansion 
in powers of m~^. 

Imposing only the Lagrange multiplier equation of mo- 
tion, the lowest order (C (to°)) piece of Eq. (58) is, 

= /3i y d^x{^F,,F,, - (d^vo)' + (d^v,)^ 
- {diVo + doVi)Fio^ 
- {f3i +p2 + Ps) J d^x (^{dovo)^ - {diVi)^) , (60) 

since V3 = 0{m~^). When f3i+P2 + p3 7^ 0, we find that 
the mode described by the dispersion relation (21) has 
energy density (at lowest order in perturbations), 

while the modes that propagate with dispersion relation 
co^ = k'j_ + have energy density, 

Ei^^/V = 2(3iiv^kl^^ and iJ^/F = 2/3ifc|fci^i, 

where ip and 62 are coordinate-independent coefficients 
of the plane wave solutions and V is the volume of space. 
Thus for (3i > the mode with dispersion relationship 
(21) has negative energy while the others have positive 
energy. The existence of a negative energy mode implies 
that the background field configuration, u'^ = mq^^ , is 
unstable when /3i -h /32 + /?3 ^ 0. 

However, as discussed in section III, when /3i + /32 + 
/33 = 0, the mode that we have just shown to have 
negative energy (when /3i > 0) vanishes. Indeed, if 
/?! + /32 + /33 = 0, then the lowest order energy functional. 




+ {d3V^)^ + {d3Vi)^), (61) 

is manifestly positive. 

VI. CONCLUDING REMARKS 

We have studied the small fluctuations about a spa- 
tially homogeneous anisotropic inflationary background. 
The anisotropy was caused by a dynamical four-vector 
that was constrained to have a spacelike invariant norm. 
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Prom the first order stability analysis of the equations 

of motion wc derived the constraints /3i > and /3i +/32 + 
/33 > for the parameters in the vector Lagrangian (1). 
Moreover, we find that all modes have positive energy in 
fiat-space and do not grow with time in the case where 
the kinetic term for the four-vector corresponds to a field 
strength tensor squared {i.e. for Pi + P2 + 03 — 0)- A 
negative energy mode propagates ii Pi + 02 + 7^ 0, 
which implies that the background given by Eqs. (2)- (5) 
is unstable when /3i + /?2 + /?3 7^ 0. 



Ae4 = Aha = -2 (Ha2Jo{t, k^) + ihj,{t, kS) , 

(A.4) 

AOs = - ^6ij^ Ahij = -ikiMt, h), (A.5) 

Aee = ^Ahis 
k± 

= {S'ik3)fo{t,k^)+if3{t,k^)) , (A.6) 



Post-submission note 



After the posting of this paper, in Ref. [23], an in- 
stability in the theory with a fixed-norm spacelike vector 
field with a field strength squared kinetic term was found 
when the physical momentum of perturbations is compa- 
rable to the Hubble scale. We considered only physical 
momenta that were much greater than or much less than 
the Hubble scale. 



APPENDIX: GAUGE ARTIFACTS 

A gauge transformation of the form in Eq. (25) induces 
the transformation k) h^i,{t, k) + Ah^i,{t, k) in 

the comoving Fourier-transformed metric perturbations. 

Exphcitly, (i, j 1, 2) 

A/ioo(t,fc) =29t/o(i,fe±), 

Ahoi{t, k) = i—^fo{t, k±) - a{t)dtfi{t, k±), 
a[t) 

Ahosit, k) = bit) (iS'ik3)dtfoit, k^) - dth{t, k^)) , 
A/iy (t, k) = -2 {Ha5,jfQ{t, k^) + ik^,fj){t, fc^)) , 

Ahiz{t, k) = -^/Ci (5'{k3)h{t, k^) + ifsit, fcx)) , 

where k± — (fci, fe) and /^(t, fcj^) are the Fourier trans- 
forms of ffj, {t, xi,X2). The corresponding transformations 
of the fields as defined in Eq. (49) are. 



and, 



Aei = Ahoo = 2dtfo{t,ki_), 



(A.1) 



Ae2 = ^Ahoi = z^/o(t,fc±) - a{t)^dtMt,k^), 

(A.2) 

Aes = A/103 = b{t) {i5'{kz)dth{t, k±) - dtfait, k±)) 

(A.3) 



A*i = '-^Ahoi 

k± 

k[iij)kkk 



A*2 



A^-a = 



kl 

€44 kq 



-a{tf-^dtfi{t,k±), 
Ahij = icijkifjit, k±), 



k± 



Ahis = 0. 



(A.7) 
(A.8) 
(A.9) 



We see that ^^(t) is invariant under the residual gauge 
transformations in Eq. (25). From Eqs. (A.2), (A.3) 
and (A.7), it is clear that one can gauge away 63 (i) 
by fixing dtf3{t,k±) and dtfo{t,k±) = 0, and one can 
gauge away Q2{t) and ^i{t) by fixing dtkifi{t,k±), and 
dteijkjfi{t,k±), respectively. 
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